Fine structure in the electronic density of states of Al-Ni-Co decagonal quasicrystal 

from ultrafast time-resolved optical reflectivity 
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We measured the temperature and fluence dependence of the time-resolved photoinduced optical 
reflectivity in a decagonal Al7i.gNin.iCoi7.o quasicrystal. We find no evidence for the relaxation 
of a hot thermalized electron gas as observed in metals. Instead, a quick diffusion of the hot 
nonthermal carries ~ 40 nm into the bulk is detected enhanced by the presence of a broad ~1 eV 
pseudogap. From the relaxation dynamics we find an evidence for the presence of a fine structure 
in the electronic density of states around ~13 meV from the Fermi energy. The structure is related 
to a weak bottleneck for the carrier relaxation observed at low temperatures. 



The ultrafast relaxation of hot photoexcited electrons 
in systems with suppressed density of electronic states 
(DOS) near the the Fermi energy (Ep) such as semimetal- 
lic graphite the cuprates[J and charge-density wave 
systems Q has received a lot of attention in the past two 
decades. The presence of a gap, or a pseudogap near Ep 
leads to a bottleneck in the hot-carriers relaxation and 
the emission of hot optical phonons.[l,[5| 

A suppression of the DOS near Ep is characteristic 
also for quasicrystals (QC).@ The pseudogap is nar- 
row in icosahedral QC and of the order of ~leV wide 
in decagonal QC (d-QC). However, there is some evi- 
dence that a finer structure might exist in the DOS near 
Ep also in d-QC. 7] In addition d-QC are characterized 
by strong anisotropies of properties along periodic axis 
and quasiperiodic (QP) plane. 0] It is therefore an in- 
teresting fundamental question what are the relaxation 
pathways of hot photoexcited electrons in d-QC which 
show metallic conductivity along the periodic direction 
and nonmetallic one in the quasiperiodic plane 0]. In 
this Letter we present results of optical pump-probe spec- 
troscopy in a decagonal AI71.9N111.1C017.0 QC. From the 
relaxation dynamics we find evidence for the presence of 
a structure in the DOS on a scale of ^13 meV around 
the Fermi energy and no evidence for the relaxation of a 
hot thermalized electron gasflfj} 

The preparation and characterization of a single grain 
decagonal AI71.9N111.1C017.0 quasicrystal was described 
elsewhere, [llj. Optical measurements were performed on 
a side facet, which is parallel to the 10-fold axis of the 
prism shaped sample mounted in an optical cryostat. A 
standard pump-probe setup was used with linearly po- 
larized pump beam with the photon energy 1.55 eV, the 
pulse length 50 fs and repetition frequency 250 kHz. The 
pump beam was focused to a 250-/im diameter spot on 
the facet in a nearly perpendicular geometry. To de- 
tect the photoinduced reflectivity AR/R a weaker probe 
beam with the photon energy, hw pro b e , either 1.55 or 3.1 
eV and the diameter 220 (im was focused to the same 
spot and detected upon reflection by a PIN photodiode. 
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FIG. 1: The probe polarization and photon energy depen- 
dence of photoinduced-reflectivity transients at 5K (a). The 
scans are normalized to enable easier comparison. Tem- 
perature dependence of photoinduced-reflectivity transients 
(b) and fluence dependence of photoinduced-reflectivity tran- 
sients at two different temperatures (c) and (d) at 1.55-eV 
probe photon energy and polarization along the periodic di- 
rection . 



In Fig. la we show the dependence of the AR/R tran- 
sients on the probe polarization at two different photon 
energies at the temperature 5K. We can identify three 
distinct timescales on which the transients show different 
behaviour. On the sub-ps timescale one can clearly iden- 
tify a sub-ps component with fast ^200-fs rise and decay 
times. The component is present in all traces except for 
the polarization along QP direction at 3.1eV-PPE where 
another component is revealed, which slowly rises on a 
few-ps timescale after initial ~100-fs rise. At an inter- 
mediate timescale between ~1 and ~20 ps we observe 
further rise of the AR/R which only starts to relax to- 
wards equilibrium after ~2 ps with 3.1-eV PPE and after 
~10 ps with 1.55- PPE. The dispersion of the AR/R rise 
indicates a two component relaxation on this timescale 
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as well. On a long timescale beyond ~20 ps all curves, 
when properly scaled, fall onto each other indicating a 
common origin of the relaxation. 

We should also note that we observed no anisotropy 
with respect to the 1.55-eV pump-beam polarization ori- 
entation on all timescales. 

To study further the temperature and pump fluence 
{T) dependence we chose 1.55-eV PPE and the polariza- 
tion along the periodic direction which has the largest 
signal to noise ratio and contains all relaxation compo- 
nents. In Fig. lb we show the temperature dependence 
of the transients at the lowest T . The sub-ps component 
is clearly temperature independent, while at the interme- 
diate timescale a shift of the AR/R peak towards shorter 
delay and decrease of the peak amplitude are observed 
with an increasing temperature. 

The F dependence at 5K (Fig. lc) is somewhat similar 
to the temperature dependence, except that at the high- 
est two T the sub-ps component shows saturation and an 
increase of the relaxation time. At 290K (Fig. Id) the 
increase of the T only influences the sub-ps component, 
while the relaxation beyond ~1 ps is ^-independent. 
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FIG. 2: Fits of equation (1) to AR/R at T = 320 ^J/cm 2 
(a) and at T — 32 /xj/cm 2 (b). For comparison a function 
proportional to 1/ %/i in (a) and a single exponential fit in (b) 
are shown by thin solid lines. The inset shows the temperature 
dependence of the fit parameter. 

We start the analysis by considering the long timescale 
first. The universal delay dependence at different polar- 
izations and PPE indicates that on this timescale all the 
subsystems are in a local thermal equilibrium and can 
be described by a single local temperature. The time 
dependence of the local temperature is governed by the 
ID heat diffusion equation. Since the optical penetration 
depth at the probe photon energies is estimated[17fl to be 
~ 25 nm, we assume that AR/R is in the lowest order 
proportional to the surface temperature rise AT s (t). 

If, after a local temperature is established, the energy 
is distributed within the surface layer of depth zq, the 



time dependence of the surface temperature is given by: 
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where D represents the heat diffusivity along QP direc- 
tion in our experimental geometry. 18J Equation (1) fits 
the measured AR/R on the long timescale rather well(see 
Fig. 2) except for low T at low T, where the temperature 
dependence of D and nonlinear dependence of AR/R on 
AT S become important. From the result of the fits we 
estimate the energy deposition depth zq and the initial 
temperature rise ATo using published values of the heat 
capacity [13, Il3( and thermal conductivity along the QP 
direction [8|, |l2j. The estimated zq shown in Fig. 3a first 
increases and then drops slowly with increasing T to ~ 
40 nm at 290K. Although we obtain at 5K virtually the 
same zq from measurements at two different T the val- 
ues below ~100K should be taken with caution due to 
the limited validity of equation (1). 

Surprisingly, at all temperatures z is significantly 
larger than the estimated optical penetration depth of 
~ 25 nm indicating a quick diffusion of hot carriers 
into the sample during the first few picoseconds of re- 
laxation. This is rather unexpected due to the strong 
Drude damping in the QP plane [Ij. One can roughly es- 
timate the hot electron diffusion constant by neglecting 
the energy transfer to the lattice, D c \ — X c \/c ip, where 
A c i is the electronic heat conductivity, c c \ the electronic 
heat capacity and p the mass density. We calculate A e i 
from resistivity [ll(. which is almost T-independent, us- 
ing Wiedemann-Franz's law. By taking c c \ from[8( we 
obtain D c \ ~ 1 cm 2 /s, which is virtually temperature 
independent. In 1 ps this gives 20 nm diffusion length 
which is of correct magnitude to explain observations. 
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FIG. 3: The temperature dependence of the initial energy- 
distribution depth (a) and the initial surface-temperature rise 
at 32-/iJ/cm 2 excitation fluence(b). 

We proceed with the analysis of the intermediate- 
timescale. Temperature dependencies of the peak am- 
plitude and the intermediate-timescale rise time shown 
in Fig. 4 clearly suggest that there exists a well defined 
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energy scale of the order of ^100K associated with the 
processes occurring on this timescale. We assume that 
AR is a consequence of the excited state absorption [3], 

Ai?cx [ deN(e)\M(e,u)\ 2 [f(e)- f T (e)], (2) 

where hw is the energy of the photons, e the energy mea- 
sured from the Fermi energy, iV(e) is the effective density 
of the initial states close to the Fermi energyflij. M(e, lu) 
the effective dipole transition matrix element, /(e) the 
non-equilibrium electron distribution function and /t(e) 
the equilibrium Fermi function. If N(e) has in addition 
to a smooth background a narrow peak at the energy E p 
there will exist a contribution Ai? p to AR proportional to 
/(-Bp) — fr(E p ) in addition to a weakly temperature de- 
pendent contribution Ai?b from the smooth background. 

In the simplest case one can assume that the nonequi- 
librium distribution is thermal, /(-Bp) = /t'(-B p ), at an 
elevated temperature which is the same as the lattice 
temperature T' = T + ATo, where ATo is estimated form 
the long timescale behaviour. In Fig. 4a we show the dif- 
ference Af T (E p ) = f T ,(E p ) — /r(-Bp) for three different 
-Bp using values of ATo from Fig. 3b. It is clearly seen 
that A/x(B p ) at E p /kB = 150K reproduces the mea- 
sured temperature dependence rather well if one assumes 
that Ai?b is temperature independent. 
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FIG. 4: Temperature dependencies of the intermediate- 
timescale peak magnitude of AR/R (a) and the intermediate- 
timescale rise time (b) at the two different fluences. In (a) 
the small open symbols connected with dotted lines represent 
A/t(-E p ) discussed in text. 

The alternative scenario is the existence of a bottleneck 
due to the suppressed density od states near the Fermi 
energy similar to the pseudogap in the cuprate supercon- 
ductors. In this case for small excitation densities and 
strong bottleneck the density of the photoexcited quasi- 
particles is given byjij], 

n po cx l/[l + Sexp(-A g /fc B T)], (3) 



where A g represents the pseudogap width and B the ratio 
between number of phonon degrees of freedom and num- 
ber of electronic states in the phonon-energy range. A fit 
of equation (3) to the data is shown in Fig. 4a. As in the 
previous case we assumed that AR is a sum of a temper- 
ature independent part Ai?b and a part proportional to 
equation (3). The resulting pseudogap width is virtually 
the same as in the previous case, Ag/fce =140K±40K. 
The large error bar comes from indeterminacy of Ai?b- 
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FIG. 5: Fluence dependence of the intermediate-timescale 
peak magnitude of AR/R at 5K. For comparison A/t(-B p ) 
is shown with E p /kB = 150K as the solid line assuming 45- 
nm initial energy deposition depth and published values of 
the heat capacity data[l, H3 , Il3l ] and reflectivity]^ in similar 
decagonal quasicrystals. 

Both models are based on an increased density of states 
at ^13 meV away from the Fermi energy and are just two 
limiting cases representing either no or a strong bottle- 
neck. In the strong bottleneck case AR pea k/R is pro- 
portional to T at small T while at large T a crossover 
to AR pea k/R oc \[T is expected. [4] We do not observe 
any of the dependencies in our case. On the other hand, 
the T dependence of AR pea k/R for small T in the no- 
bottleneck model is compatible with the experimental 
one shown in Fig. 5. Unfortunately, the no-bottleneck 
model fails to fit the data beyond T « 1 mJ / cm 2 as well. 
The reason for this might be in ^"-dependence of Ai?b, 
which also shows some saturation beyond T ~ 3 mJ / cm 2 
judging from the measurements at 290K.(see Fig. Id). 

While the JF-dependence clearly excludes the strong 
bottleneck scenario it does not exclude a weak bottleneck, 
the presence of which is suggested by the increase of the 
intermediate-scale risetime below 100K (see Fig 4b). In 
order to elaborate this we need to analyze the response 
on the shortest timescale first. 

The sub-ps peak amplitude has linear ^-dependence 
and the decay time has no ^-dependence up to T f=a 3 
mJ/cm 2 . In addition, both show no T-dependence so 
the sub-ps dynamics cannot be attributed to the thermal 
relaxation [l0( where /(e) is assumed to be the Fermi dis- 
tribution with the nonequilibrium electronic temperature 
T e . /(e) is therefore nonthermal on the sub-ps timescale 
and the sub-ps peak is due to the relaxation across the 
broad pseudogap [rj [l5| and the hot carriers diffusion into 
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the sample. The broad pseudo-gap suppresses the prob- 
ability for the creation of the low energy electron-hole 
(e-h) pairs increasing the electron thermalization time 
beyond ~500 fs enhancing the hot carrier diffusion into 
the sample. 

There is also no evidence for the thermal relaxation of 
the hot electron gas on a timescale up to several ps. The 
intermediate-timescale risetime increases with decreas- 
ing temperature which, is inconsistent with the thermal 
reraxationa where the relaxation time is predicted to 
be proportional to the lattice temperature Tl for T c ~ 
T"l- The intermediate-timescale dynamics might there- 
fore at least in part be attributed to the dynamic lattice 
expansion due to an increasing lattice temperature fljj]. 
The expansion, however, also cannot explain the increase 
of the intermediate-timescale risetime at low tempera- 
tures, implying another relatively slow relaxation chan- 
nel. We tentatively assign this channel to the relaxation 
of hot nonthermal optical phonons. We believe that due 
to the increased thermalization time a significant num- 
ber of hot optical phonons is generated in addition to the 
low energy e-h pairs during the thermalization. These 
nonthermal hot optical phonons can decay by the an- 
harmonic decay or by excitation of the e-h pairs. At 
low temperatures the anharmonic decay channel is sup- 
pressed resulting in a bottleneck and the increase of the 
intermediate-timescale risetime below 100K. 

The relaxation upon photoexcitation in 
AI71.9N111.1C017.0 quasicrystal therefore proceeds 
roughly in two steps. In the first step the presence of 
the broad pseudogap slows down the hot carrier ther- 
malization. During the relatively slow thermalization 
hot carriers diffuse up to ~40 nm into the bulk and 
a part of the absorbed energy goes to the nonthermal 
hot phonons. During the second step, which is at 5K 
characterized by the relaxation time of ^3 ps, the hot 
nonthermal phonons decay. At low temperatures the 
nonthcrmal-phonon decay is suppressed due to the 
decrease of the anharmonic relaxation rate and the 
weak bottleneck in the relaxation through the electronic 
channel. The bottleneck is a consequence of a fine 
structure in the electronic DOS at ~13 meV from the 
Fermi energy. 
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